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In the framework of a two-band model, we study the phase separation regime of different kinds of 
strongly correlated charge carriers as a function of the energy splitting between the two sets of bands. 
The narrow (wide) band simulates the more localized (more delocalized) type of charge carriers. By 
assuming that the internal chemical pressure on the Cu02 layer due to interlayer mismatch controls 
the energy splitting between the two sets of states, the theoretical predictions are able to reproduce 
the regime of phase separation at doping higher than 1/8 in the experimental pressure-doping- Tc 
phase diagram of cuprates at large microstrain as it appears in overoxygenated La2Cu04. 

PACS numbers: 74.72.-h,71.27.+a, 64.75.-g 



I. INTRODUCTION 

The mechanism driving the emergence of a quantum 
macroscopic coherent phase that is able to resist to the 
de-coherence effects of high temperature remains a major 
topic of research in condensed matter. The realization 
of this macroscopic quantum phase in doped cuprates 
close to the Mott insulator regime has stimulated a large 
amount of investigations on the physics of strongly cor- 
related metals. Most of theoretical papers treated mod- 
els of a homogeneous system made of a single electronic 
band (or models of multiple hybridized bands reduced to 
a single effective band), with a large Hubbard repulsion. 

There is growing agreement that the solution of the 
problem of high- Tc superconductivity requires the correct 
description of the normal state where spin, charge, orbital 
and lattice degree of freedoms compete and the functional 
phase emerges in a complex system with two main com- 
ponents showing mesoscopic phase separation. Here we 
consider a theoretical model of the mesoscopic phase sep- 
aration in a two-band scenario of two strongly correlated 
electronic fluids. This simple model grabs the key physics 
of the anomalous normal phase in cuprates exhibiting the 
phase separation as a function of charge density and the 
energy splitting between the two bands. This allows the 
understanding of the different superconducting phases in 
different cuprate families, i.e., the new 3D phase diagram 
where the critical temperature depends on the doping 
and pressure^. The motivation of this theoretical work 
is based on the results of recent experiments using angu- 
lar resolved photoemission spectroscopy (ARPES)^i2ii^'^ 
and scanning tunnelling spectroscopy (STM)ii^ provid- 
ing compelling experimental evidence for dual nature of 



charge carriers and the nanoscale phase separation of 
the two components in two different spatial domains in 
cuprate high-Tc superconductora^ i^^'^^ . 

A clear case for the phase separation of the two types 
of charge carriers is overoxygenated La2Cu04+j, where 
the interstitial oxygen ions are mobile above 180 K and 
stimulate the phase separation of the two different kinds 
of dopant hole9i^>i^>Hii5a6a7j^j^^^ Currently, from 

the analysis of magnetic neutron scattering experiments, 
there is an agreement for the frustrated mesoscopic phase 
separation at doping larger than 1/8 in Sr-doped La214, 
Y123, and Bi2212 between a first more delocalized com- 
ponent that does not show spin fluctuations and a sec- 
ond more localized electronic component, showing stripe- 
like spin fluctuations2ii22. Several reviews and books 
have been published on the two-component scenario and 
phase separation in cuprate a-'^^'^^i^^'^'^i^^'^^i^^ . Here we 
focus not on the well studied phase separation in the 
underdoped regime, near the Mott phase, between a 
hole-poor antiferromagnetic phase and a metallic hole- 
rich phase but on the phase separation in the over- 
doped regime^^^^ii^i^^^^ between a hole 
poor phase with doping close to 1 /8 and a hole-rich phase 
with doping close to 1/4. 

This phase-separation scenario has been described 
in Refs. [T8ll2lll2^ for LaSrCuO and overoxygenated 
La2Cu04. It is based on the experimental fact that vari- 
ation of the magnetic incommensurability due to s pin 
stripes saturates at doping 1/8 see Fig. 15 of Ref. [2ll . 
The residual magnetic scattering at high doping suggests 
that one of the phases has stripe correlations similar to 
the a; = 1/8 phase, with the volume fraction of t his p hase 
decreasing with indicated in Fig. 17 of Ref.l2ll. The 
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other phase is presumably uniformly doped. The picture, 
then, is that as one increases the doping beyond x = 1/8, 
it becomes unfavorable to accommodate the additional 
holes in stripes; instead, patches of the uniformly-doped 
phase grow at the expense of the stripe phase. The max- 
imum Tc seems to occur in a mixed phase region domi- 
nated by the stripe phase. 

A similar scenario is now well accepted for understand- 
ing the physics of phase separation in manganites22i2^i2^. 
It was shown that even in the absence of any specific 
order parameter, the presence of two strongly corre- 
lated electron bands leads to the possibility of a phase- 
separated staters. 

In Ref. [13, the evolution of phase separation was stud- 
ied as a function of doping. However, a large amount of 
data clearly indicates that the phase separation regime is 
not only a function of doping but also of the anisotropic 
chemical pressure acting on the Cu02 layers, due to in- 
terlayer mismatch^i^^i^i^^. The chemical pressure is a 
well established physical variable that controls the phys- 
ical properties of perovskites and it is usually measured 
by the average ionic radius of the cations in the inter- 
calated layers or the tolerance factor t, in fact, the in- 
ternal chemical pressure in perovskites can be defined as 
T] = 1 — t. In all perovkites and particularly in mangan- 
ites, it is well established that the phase diagram of the 
electronic phases depends on the two variables, charge 
density and chemical pressureS^. Since the early years of 
high-Tc superconductivity research the mismatch chem- 
ical pressure has been considered as a key variable con- 
trolling the electronic properties of cuprates only on one 
family, La214'^^, however it was not possible to extend 
this idea to other families for the presence of a plurality 
of intercalated layers with cations having largely differ- 
ent coordination numbers. Therefore, it was not possible 
to compare the average ionic size (va) in the intercalated 
layers and to get the tolerance factor t for all cuprate 
families . This problem was solved by obtaining the in- 
ternal chemical pressure from the measure of the com- 
pressive microstrain e — (Rq — r)/r in the Cu02 plane 
(that has the same absolute value as the tensile micros- 
train in the intercalated layers) where r is the average 
Cu-0 distance and i?o — 0.197 nm is the unrelaxed Cu- 
O distance^i^i^. Therefore, the chemical pressure is 
proportional to microstrain, rj — 2e. 

In this new 3D phase diagram, the phase separation for 
the overdoped regime in overoxygenated La214 occurs in 
a family with high chemical pressure close to rj — 8%, 
while it becomes a frustrated phase separation in the 
LaSrCuO, Bi2212, and Y123 that are in the range of 
chemical pressure 7% > rj > 4%, while for cuprates with 
lower microstrain only very fast critical fluctuations could 
be present^i^'^. In this paper, we propose a model of 
a two-component system made of two different strongly 
correlated electron bands where the chemical pressure 
controls the energy splitting between the two bands. The 
phase separation in the overdoped regime can exist for 
specific values of the ratio between the bandwith of the 



two bands. The critical point for the transition from 
a frustrated phase separation to a non frustrated phase 
separation can be obtained by tuning the long range 1/r 
Coulomb repulsion that frustrates the phase separation 
as going from Sr doped to oxygen doped Lal24. 



II. THE MODEL 

The existence of the two types of the strongly cor- 
related charge carriers in cuprates can be described in 
terms of the two-band Hubbard model. The Hamilto- 
nian of such a system can be written as^ 

{nm)a,a- ncr na,cr 

1 U' 

^ ^ ^ ^ U l^naa'^naa ^ ^ ^nacr^nacr' ■ (1) 

na,(T nOL.aa' 

Here, a\^^„ and the creation and annihilation op- 

erators for electrons corresponding to bands a = {a, 6} at 
site n with spin projection tr, and n-aaa = <^\iaa^-naa- The 
symbol (...) denotes the summation over the nearest- 
neighbor sites. The first term in the right-hand side of 
Eq. Ilj corresponds to the kinetic energy of the conduc- 
tion electrons in bands a and h with the hopping integrals 
ta > tfc. In our model, we ignore the interband hopping. 
The second term describes the shift Ai? of the center of 
band h with respect to the center of band a (A_E > if 
the center of band h is below the center of band a). The 
last two terms describe the on-site Coulomb repulsion of 
two electrons either in the same state (with the Coulomb 
energy J7") or in the different states (J7'). The bar above 
a ox a denotes not a or not a, respectively. The assump- 
tion of the strong electron correlations means that the 
Coulomb interaction is large, that is, U", U' ^ t^, AE. 
The total number n of electrons per site is a sum of elec- 
trons in the a and b states, n = Ua + ni,, and fi is the 
chemical potential. Below, we consider the case n < 1 
relevant to cuprates. 

Model H]) predicts a tendency to the phase separation 
in a certain range of parameters, in particular, in the 
case when the hopping integrals for a and b bands differ 
significantly {ta > t^)^^. This tendency results from the 
effect of strong correlations giving rise to dependence of 
the width of one band on the filling of another band. In 
the absence of the electron correlations (n ^ 1), the half- 
width Wa = zta of a band is larger than Wb = zti, {z is 
the number of the nearest neighbors of the copper ion). 
Due to the electron correlations, the relative width of a 
and b bands can vary significantly'^'^ . 

In the limit of strong correlations, [/",[/' — > oo, we 
can describe the evolution of the band structure with the 
change of n and A_B following the method presented in 
Ref. [so. We introduce one-particle Green function 

Ga,(n - no, t " to) = -i(f an„,(i)a|,„„,(io)), (2) 
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where T is the time-ordering operator. The equations of 
motion for one-particle Green function with Hamiltonian 
(H]) includes two-particle Green functions 

^acr,/3cr'(n-no, t - t^) = -i{fanaa{t)nnPa' {t)a'l^_^^„{tQ)) . 

In the considered limit of strong on-site Coulomb re- 
pulsion, the presence of two electrons at the same site is 
unfavorable, and the two-particle Green function is of the 
order of 1/U, where U ^ Ua, U' . The equation of motion 
for Gaa.fJa' includcs the three-particle terms coming from 
the commutator of a^aait) with the U terms of Hamil- 
tonian ([1]), which are of the order of and so on. 
In these equations, following the Hubbard I approach"^-, 
we neglect the terms of the order of and make the 
following replacement (Tan+maa(t)n„/3<T' (O«noaa(*o)) ^ 



("-n/3<T')(ra„^ 



r(to)). As a result, we derive 
a closed system for the one- and two-particle Green 
functions2^i2i,. This system can be solved in a conven- 
tional manner by passing from the time-space {t, r) to the 
frequency- momentum (w,k) representation. In the case 
of superconducting cuprates the total number of electrons 
per site does not exceed unity, n < 1. The upper Hub- 
bard sub-bands are empty, and we can proceed to the 
limit Ua, U' oo. In this case, the one-particle Green 
function Gaa is independent of U and can be written in 
the frequency-momentum representation asS2i^ 



LU + 1-1 + Ai?" - gacrWaCO^) ' 

where AE"' = for a = a and A^;" = AE for a = b 



(3) 



9a 



1 



(4) 



and po{E') = J d^kS{E' ~ C(k))/(27r)3 is the density of 

states for free electrons (with the energy normalized by 
unity, \E\ < 1). The chemical potential /i in Eq. ([5]) can 
be found from the equality n — Ua + rib. 



Wulw = 
'"^b a 
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Figure 1: (Color online) Evolution of the occupation numbers 
Ua and Wi, of the bands a and 6 at fixed doping 5 = 1 — 
n = 0.3 in the absence of phase separation. The region of 
phase separation lies between two vertical dotted lines. There 
we have two phases: Pa, including mostly a charge carriers 
and Pf, with dominant b carriers. The content of diflterent 
types of carriers in Pa and Ph is given by the intersections of 
Ua and Uh curves with left and right dashed vertical curves, 
respectively. The change in concentration p of phase Pa in the 
phase-separation region is shown by the (green) dot-dashed 
line. 



naa = {nnaa) IS the average number of electron per site 
in the state {a, a), and C(k) is the spectral function de- 
pending on the lattice symmetry. Since the results does 
not vary crucially with the change of the lattice symme- 
try^, here we consider the case of the simple cubic lat- 
tice, when C(k) = - [cos(A;^d) -I- cos(fc^d) -I- cos(fc^d)] /3, 
d is the lattice constant. In the main approximation in 
1/U, the magnetic ordering does not appear and we can 
assume that — = na/2. 

Equations ([3]) and (|4]) demonstrate that the filling of 
band a depends on the filling of band b and vice versa. 
Indeed, using the expression for the density of states 
p^(E) = -TT-^Im / Ga{k,E + iO)d^k/{2TT)^, we get the 
following expression for the numbers of electrons in bands 
a and b (a = a, b) 



2gano 



p + AE°' 

gaWa 



where 



no{fi')= / dE'poiE'), 



(5) 



(6) 



When the energy band b is far above the center of band 
a {AE < 0), there exist only a electrons. With the in- 
crease of AE, the chemical potential reaches the bottom 
of the b band —AE — Wb- At higher AE, the b electrons 
appear in the system, and the effective width of a band, 
= 2waga, starts to decrease. At large positive values 
of AE, the a carriers in the system disappear and there 
exist only b electrons. The plots of ria, rib, and the effec- 
tive bandwidth as functions of AE are shown in Figs. [1] 
and [21 respectively. 

The energy of the system in homogeneous state, i^hom, 
is the sum of electron energies in all filled bands. Simi- 
larly to Eq. ([5]) , we can write i^hom in the form 



AE"^ 



gaWa 



where 



dE'E' po{E') . 



AEub, (7) 



(8) 



The dependence of -Ehom(^) is shown in Fig. [3] at dif- 
ferent values of AE. We see that within a certain n 
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Figure 2: (Color online) Effective widths w^fj of the a and 
b bands versus band shift AE at fixed doping 5 = 1 — n = 
0.3. The (green) dot-dashed curve illustrates the behavior of 
the chemical potential /x; the (green) hatched area under this 
curve corresponds to the states occupied by charge carriers. 




Figure 3: (Color online) The energy of the system vs doping 
level n at different values of AE. Solid curves correspond to 
the homogeneous state, whereas the dashed curves are the 
energies of the phase-separated state without taking into ac- 
count electrostatic and surface contributions to the total en- 
ergy. 



the total energy due to surface effects and a charge redis- 
tribution. However, at first we do not take into account 
this effects. We consider two phases. Pa (low carrier den- 
sity) and Pb (high carrier density), with the number of 
electrons per site ni and rt2, respectively. A fraction p of 
the system volume is occupied by the phase Pa and I —p 
is a fraction of the phase Pb- We seek a minimum of the 
system energy 

E°^{ni,n2) = P-Ehom(«l) + (1 - p)-Ehom('^2) (9) 

under the condition of the charge carrier conservation 
n = pni -|- (1 — p)n2. The results of calculations of the 
system energy in the phase-separated state are shown 
in Fig. [3] by the dashed lines. We see that the phase 
separation exists in the range of n where both types of 
charge carriers coexist in the homogeneous state. The 
ratio of the numbers of a and b carriers is different in 
different phases. In the first phase Pa, almost all charge 
carriers are in the band a, while in the second phase Pb 
the situation is opposite. 

The redistribution of charge carriers in the phase- 
separated state gives rise to the additional electro- 
static contribution, Ec, to the total energy. This term 
in the Wigner-Seitz approximation was calculated in 
Ref. [10. At p < 0.5 it can be written as Eq = 
V (rii — 712)^ {Rs/d)'^ u{p), where 



u{p) ^ 2Trp (2 - 3p^^^ + /5, 



(10) 



and V is the characteristic energy of the intersite 
Coulomb interaction and Rg is the radius of the spher- 
ical droplet of the phase Pa surrounded by the shell of 
the phase Pb- In the case p > 0.5, we should replace 
ni <-!■ 712 and p <-> 1 — p. The second contribution to the 
total energy, depending of the size of inhomogeneities, 
is related to the surface between two phases. The cor- 
responding energy per unit volume can be presented in 
the form Eg — pSa{ni,n2)/VQ, where p < 0.5, S is the 
surface and Vq is the volume of the inhomogeneity, and 
(j(ni, 712) is the surface tension. For e spherical droplets, 
we have Eg = 3pa{ni,n2)d/ Rs- lip > 0.5, we should re- 
place p — > 1— p. Minimization of the sum Eqs = Eq + Eq 
with respect to Rs allows us to calculate this value. In 
doing so, we get at p < 0.5 



R.^d 



3pa(ni,n2) 
2Vb(772 - ni)'^u{p) 



1/3 



(11) 



range, the system can have a negative compressibility, 
d'^E\^ora/dn^ < 0, which means a possibihty for the 
charge carriers to form two phases with different elec- 
tron concentrations22iH. The negative compressibility 
disappears when the centers of the bands are far apart 
from each other (see, e.g., the curve corresponding to 
AE/wa = 0.4 in Fig. [3]). 

The phase separation may be hindered by increase of 



The total energy of the inhomogeneous state then reads 

Eps = pi^hom(71l) + (1 - p)£;hom(772) + EcsiRs) , (12) 

where Rg is given by Eq. The surface energy 

comes from the size quantization and it was estimated 
in Ref. [13. The electrostatic and surface contributions to 
the energy related to an inhomogeneous charge distribu- 
tion reduce the range of n, in which the phase separation 
is favorable, see the phase diagram in Fig. [D 
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III. RESULTS AND DISCUSSION 

The undoped state of the cuprates corresponds to 
one electron per site {n = 1) in the model used in 
Ref. [lO- The number of itinerant holes 6 is related to 
n as S = 1 — n. In general, the relationship between 
n and 5 could be more complicated^, however, for the 
present considerations such corrections are not of prin- 
cipal importance. The phase diagram of the model (P) 
in the {S, AE) plane is drawn in Fig. U) In this figure, 
below the lower (red) solid line, we have the charge carri- 
ers only of a type, whereas above the upper (blue) curve, 
there are only b carriers. If we ignore the possibility 
of the phase separation, the relative number of a and b 
charge carriers varies gradually between these two lines. 
The evolution of the occupation numbers and rif, of 
the two bands with AE at a fixed doping is illustrated 
in Fig. [T] Note also that the effective widths of the a 
and b bands also vary with the band shift due to the 
electron correlation effects. This is illustrated in Fig. [21 
Taking into account possible phase-separated states re- 
sults in a significant modification of the phase diagram in 
the range of intermediate doping. In the hatched (green) 
region in Fig. 3] the homogeneous state becomes unfa- 
vorable and the system separates into two phases {Pa 
and Pf,) with different numbers of charge carriers per site 
m ~ Ua and n2 ~ ni,. The electrostatic contribution 
to the energy related to an inhomogeneous charge dis- 
tribution reduces the doping range, in which the phase 
separation is favorable^S. In Fig. [H we illustrate that a 
relatively small energy loss due to the charge dispropor- 
tionalization leads to a substantial decrease in the area 
of the phase-separation region: compare the areas indi- 
cated by arrows corresponding to V/wa = 0.01 and 0.015 
{V is the characteristic energy of intersite Coulomb in- 
teraction'^*) and the whole hatched area corresponding to 
V ^ 0. Note that at low hole doping (n close to one), the 
antiferromagnetic (AF) correlations are dominant, which 
requires a special analysis. 

Estimating the contribution of the long-range Coulomb 
interaction to the energy of the phase-separated state, 
we assume the simplest droplet-like geometry of the in- 
homogeneities. The phase separation occurs in the range 
of parameters where the energy of homogeneous state 
as a function of doping has a negative curvature cor- 
responding to the negative compressibility2Si2£. It was 
widely discussed in literature that long-range Coulomb 
interaction in the systems with negative compressibility 
can give rise to more complicated geometry of the phase 
separation (stripes, layers, rods, etc.), see Refs. ISQlliol 
and references therein. Thus, a due account of the 
long-range Coulomb interaction could reproduce differ- 
ent superstructures (stripes, in particular) observed in 
the cuprate superconductors near the optimum doping. 
However, the proper analysis of the inhomogeneity geom- 
etry requires a further study based on a more complicated 
model. 

Now let us discuss the relation of the above model 
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Figure 4: (Color online) The phase diagram of model ((T} in 
the (doping, band shift) plane at the ratio of band widths 
Wb/wa = 0.3. Below the lower (red) solid line, there are 
charge carriers only of a type, whereas above the upper (blue) 
curve - only b carriers. Between these lines, there appears the 
region of phase separation (marked by (green) hatching) . The 
charge disproportionalization in the phase-separated state can 
substantially reduce this region: the arrows indicate the phase 
separation regions at V/wa = 0.01 and 0.015, where V is the 
characteristic energy of intersite Coulomb interaction. 

to the experimental situation in the copper-based per- 
ovskites, where two types of charge carriers and inho- 
mogeneous (phase-separated) state are observed. The 
inhomogeneous state in cuprates corresponds to the co- 
existence of two phases. One of them is characterized by 
a superstructure (charge ordering, stripes, etc.) and an- 
other one has no superstructure. The state with charge 
(or spin) superstructure corresponds to a higher degree of 
localization and, therefore, to a smaller value of the hop- 
ping integral. Naturally, a charge carrier may hop either 
retaining short-range order and gaining in the potential 
energy or hop in an arbitrary way with larger hopping 
integral thus gaining in the kinetic energy. The former 
corresponds to our b state and the latter to the a state. 
In our analysis, we did not consider any ordering, which 
arises in the next-order approximations. In particular, 
magnetic order requires taking into account the terms 
of the order oftl^/U and the charge ordering implies al- 
lowing for the Coulomb interaction of carriers at different 
sites (nearest-neighbor at least). 

The relative position of the two bands, AE, and hop- 
ping integrals, ta.b^ depends, in particular, on the chem- 
ical pressure proportional to a microstrain e in the crys- 
tal lattice. To describe the experimental phase diagram 
of cuprates in the ((5, e) plane^ i^^i'^^'^'^ , we should know 
the relationship between the model parameters and the 
chemical pressure. It is natural to assume that the two 
bands in the cuprate crystal originate from the double 
degenerate Cg hole level of Cu-^+ (configuration (f) in 
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the crystal field of cubic symmetry. The splitting of this 
level occurs due to lattice distortions related to the Jahn- 
Teller effect, lowering the cubic symmetry. The chemical 
pressure distorts the crystal lattice even more and should 
affect the value of AE significantly. It is natural to as- 
sume that AE and e are linearly related to each other, if 
\e\ <^ 1. So, we can write 

AE{e) = AE{0) + AEif{e), (13) 

where /(e) is a dimensionless function and /(e) « e 
at |e| <C 1. Cu^+ is a typical Jahn- Teller ion and we 
can assume that AEi is of the order of the characteris- 
tic Jahn- Teller energy, which is larger than ta (see, e.g., 
Refs. ,30>38,41i and references therein). The effect of the 
microstrain on the relative band positions can be signifi- 
cant since the value AE arises due to splitting of the orig- 
inally degenerate levels. In the same time, small strains 
give rise only to a small corrections to the bandwidth. So, 
the ratio tb /ta is considered further on as independent of 
e. Note also that the values of the intersite Coulomb in- 
teraction V characteristic of perovskites is of the order 
of 0.1 — O.Olia (see, e.g., Ref . ISSl and references therein). 
Thus, the value V = O.Olbzta used below is quite reason- 
able. 

Bearing this in mind, we compare the theoretical phase 
diagram in Fig. d] with the experimental 3D phase dia- 
gram of cuprates^'^^ in Fig. O The left-hand y scale is 
the mismatch chemical pressure rj related to microstrain 
as ?7 = 2e and the x axis is the doping (the number of 
holes per Cu site) . The color plot represents the values of 
critical temperature in different superconducting cuprate 
families. The plot shows the fit of the experimental data 
of a large number of materials with the convolution of 
a parabolic curve with the maximum at Tmax for Tc as 
function of doping, and an asymmetric Lorentzian for 
Tmax as a function of the mismatch chemical pressure 
with the maximum of 135 K at 2e = 4%. The y axis in 
the right-hand side of the figure gives the energy distance 
AE between the center of band a and band b normalized 
to the width of band Wa of the more itinerant carriers. 
The phase diagram involving the superconducting criti- 
cal temperature, chemical pressure, and doping reaches 
the Tc maximum at 2e = 4% and 0.16 holes per Cu sites. 
Based on the aforementioned consideration, we can take 
AE{0) = -0.133u;a and AE{0) = 6.67wa in Eq. P^. 
We can identify a low-doping insulating phase, for any 
chemical pressure, at doping smaller than 0.06, where the 
vertical dashed line indicates the line of metal-insulator 
transition. The experimental investigations of the novel 
3D phase diagram of cuprates indicate that the homoge- 
neous metallic phase, with more delocalized states, oc- 
curs for both high doping and low chemical pressure i.e., 
in the low right corner of the figure. In this region, we 
have in the theoretical model the charge carriers only 
in the band a. On the contrary, the homogeneous phase 
made of localized states where the striped phase appears, 
occurs at the corner on the top-left side of the figure. 
In the theoretical model, we have in such a region, the 



charge carriers only in the band b. The superconducting 
phase occurs in the intermediate region between these 
two limiting cases. The phase separation region pre- 
dicted by our model (inside the area bounded by the solid 
white line) corresponds to the superconducting phase. It 
is in the qualitative agreement with the STM, extended 
X-ray absorption fine structure (EXAFS), neutron pair 
distribution functions (PDF) experiments showing that 
high- Tc superconductivity occurs in a regime of nieso- 
scopic phase separation. The present results show that 
the maximum critical temperature occurs where the en- 
ergy splitting between the more itinerant, band a and the 
more localized, band b is close to zero. 



T (K) 




Number of holes per Cu site 



Figure 5: (Color online) The values of superconducting tran- 
sition temperature Tc from for dark blue to 135 K for dark 
red is shown in a color plot as function of chemical pressure 
{2e) and doping (holes number per Cu site). The white curve 
corresponds to the phase separation region given by the two- 
band Hubbard model ([1} corresponding to wi,/wa = 0.3 and 
V/wa ~ 0.01 (see Fig. |4}. Phases Pa and Pb include mostly 
the carriers of a and 6 types, respectively. The black solid line 
is the boundary of the phase-separated state deduced from 
neutron scattering and anomalous diffraction experiments for 
cuprates (La214, Bi2212, and Y123 systems). 

The theoretical phase diagram reproduces qualita- 
tively the experimental results on the phase separation in 
cuprates (La214, Bi2212, and Y123 systems) near opti- 
mum doping obtained by neutron scattering and anoma- 
lous diffraction techniques^i^i^i. The phase separation 
arises in the intermediate doping range and disappears 
at low and high doping levels. The phase with "more 
itinerant" electrons exists at small microstrains, "more 
localized" (and more ordered) phase arises at higher mi- 
crostrains, and the phase-separated state is located in 
the intermediate range of e. However, our calculations 
predict the phase separation in a broader doping range 
than in the experiments. It seems to be a consequence of 
simplifications used in the formulation and approximate 
analysis of the two-band Hubbard model ([1]). To improve 
the agreement with the experiment it is necessary to take 
into account specific features of the lattice and electron 
structure of the cuprate superconductors. In particular, 
we disregard the interband electron transitions, that is, 
we neglect the terms tabcinaa-'^mbcr in Hamiltonian ([1]) as- 
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suming that tab = 0. The doping range where the phase 
separation can exist reduces with the increase of iair^- 

As it foUows from Fig. [H the superconducting tran- 
sition temperature Tc is the highest for the parameters 
range where the system is in the phase-separated state. 
This is an indication that the mechanism of the phase 
separation is intimately related to the phenomenon of 
the superconductivity. It is worth to note that in the case 
when interband coupling tab is in the range ta < tab < 
the electron density of states has a peak near the Fermi 
level in the parameter range corresponding to the phase- 
separated state, where Tc is maximun*^. We can not 
claim whether this fact is accidental or not. 



IV. CONCLUSIONS 

Up to now, most of the attention both of experimen- 
talists and theorists has been addressed to the phase 
separation in the underdoped regime between a first 
undoped antiferromagnetic phase and a second doped 
metallic phase of cuprates. Now, we have an evidence 
for mesoscopic phase separation in the overdoped region 
of cuprate superconductors where a striped phase at dop- 
ing 1/8 coexist with a metallic phase with doping close 



to 1/4. In our paper, we were dealing just with this sit- 
uation. 

We have presented an emerging theoretical scenario 
that relating the phase separation to the chemical pres- 
sure. This scenario grabs key physical aspects of the 
3D phase diagram of cuprates. It was shown that the 
two-band model is appropriate for the normal phase of 
all cuprate superconducting families, where the energy 
splitting between the two bands is controlled by mis- 
match chemical pressure. In the regime where the two 
bands are close in energy, the system is unstable toward 
the phase separation. The highest critical temperature 
of the superconducting transition in cuprates is attained 
within the phase-separated state. 
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